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Abstract 

We estimate in this paper the errors in the zenith angle distribution for the charged current events of 
the solar neutrinos caused by the uncertainty of the earth electron density. In the model of PREM with 
a 5% uncertainty in the earth electron density we numerically calculate the corrections to the correlation 
between [IV] 5 /[IV] 2 and [IV] 2 /[IV] 3, and find the errors notable. 

Forthcoming results from SNOQ include a measurement of the day-night asymmetry g, g. 

This measurement is crucial to confirming the matter conversion solution to the solar neutrino problem. And 
the analysis on the zenith angle distribution of the events during the night may provide some insights to 
distinguish the various MSW solutions, i.e. LMA, LOW and SMAjfJ . In the calculation of the regenerated 
v e flux, the electron density of the earth matter with which the neutrinos interact is a critical quantity. The 
uncertainty in Earth-matter density and chemical component can be a major cause of the error in Ajjn 
and the zenith angle distribution. So it will be interesting to estimate these errors. Furthermore since the 
experimental value of the Adn is around ~ 0. 047 [g and the theoretical expectations on the zenith angle 
distributions are small in magnitude^, it is necessary to perform a quantitative estimation on these errors. 

In this paper, we follow the procedure outlined in [7] and study the uncertainty in the earth matter 
density, then investigate its implications on the predictions of Aon and the zenith angle distributions. We 
quantify the uncertainties of the earth matter in terms of two parameters: one is SN e /N e , the variation 
in magnitude of the density which generally is expected to be around a few percent; the second one is 5x 
which specifies the limitation on the spacial dimension by geophysics experiments and inverting calculations 
used in the fit of the earth density models. In general the scale Sx is not much larger than the neutrino 
oscillation length, e.g., in the case with the parameters of the favored LMA solution, so its effect might arise 
beyond the linear order. We will show in this paper that this effect causes an sizable error in the zenith 
angle distributions. 

To begin with we consider a two-neutrino mixing model for simplicity. As discussed ing, g] the neutrino 
can be treated as a incoherent mixture of two mass eigenstates. In the day-time the survival probability for 
v e is given by, 

P D = P lC os 2 d + (1 - Pi) sin 2 (9 , (1) 

where the mixing angle is defined through, 


v\ = cos 6 v e — sin 6 , ^2 = sin 9 v e + cos 9 zz M , (2) 

and P\ is the probability of the v e —> v\ conversion inside the sung, g. During the night time, the presence 
of the earth matter leads to a zenith angle dependent regeneration of the u e , 


P N = P 1 + ( 1 - 2 P 1 )P 2e =P D - 2 Xf reg , 

where P 2 e is the probability of the V 2 —> v e conversion inside the Earth, X = P\ ~ 1/2. And 

freg(9z ) = P 2 e(earth matter ) — P 2 e {vacuum ) , 

is the regeneration factor which vanishes in the absence of the Earth matter effect. Defining f reg 
regeneration factor integrated over the zenith angle, one has the day-night asymmetry, 
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The matter effects has entered the day-night asymmetry through f reg - Formally, 

pD cos 6 Z 

P 2 e (E V ,0 Z )= cos 6 [T exp [-i 
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where D = 12742 is the diameter of the Earth in unit of kilometer and H[N^ Z (a;)] is the effective Hamiltonian 
for the given trajectory with zenith angle 9 Z , 
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In Eq. (0), Nl- {x) is the Earth Electron Density (EeD) along the trajectory of the zenith angle 9 Z . If the 
density is known, the regeneration factor f reg = P 2 e — sin 2 6 can be calculated accurately. As an example 
we take the Preliminary Reference Earth Model (PREM) ]To| and plot in Fig.l.(a) the regeneration factor as 
a function of the zenith angle. In the numerical calculation we take the neutrino energy to be 11 MeV and 
the oscillation parameters to bcfn| 

LMA : Ami 2 = 3.7 x 10" 5 , tan 2 6 = 0.37 , 

LOW : Amis = 1.0 x 10“ 7 ,tan 2 6 = 0.67 (8) 


One can see from this figure that the regeneration factors oscillate periodically with certain lengths. And 
different oscillation lengths correspond to different MSW solutions. 

Given the parameters in Eq.(||) and the standard solar density Q, we follow || and obtain numerically 
that cos 26s ss —1 and P c ~ 0, which can be used to get the v e —> v\ conversion in the Sun. Fluctuations 
in the solar density will affect Pi, consequently also influence the MSW solutions[p~3|]. In that situation a 
variance of P D has been defined to estimate the relevant error Jm|. In this paper, however, we concentrate 
on the errors caused by the uncertainty of EeD. 

The EeD available today, is known only to some certain precision [|1| |l(|. As to PREM, significant 
uncertainties due to the local variation have been documented ||l 7| . Quantitatively its precision is roughly 
5% averaged per spherical shell with thickness of 100 Km or so [|18| . The uncertainties of the Earth matter 
density cause errors in the calculation of the v e survival probability during the night time. In the following 
we study numerically the uncertainties in the solar neutrino zenith angle distributions. 

As described detailly in j?J , we introduce a weighted average over the whole sample space of possible earth 
density profile. Denoting the averaged earth density function, such as the widely used PREM by N e (x), we 
have N e {x) =< N e (x) >= J[DN e \F[N e (x),x][VN e \ where F[N e (x), x\[UN e \ is the probability of obtaining 
the EeD N e (x) in the neighborhood of x, 

F[N e (x),x\ = exp{ — In 2 [ N e (x)/N 0 (x) J/[2s 2 (a;)l}, 

-/V e (x)v 27 rs(x) 

s(x) = i/lnfl +r 2 (x) 1, N 0 (x) = N e (x)exp[~s 2 {x)/2 J, (9) 

where r(x) = a(x)/N e (x) characterizes the precision of the earth electron density. 

The averaged value and the variance of the zz 2 —*■ v e conversion probability can be written now separately 
as, 


<P 2 e{E v , 6 z )> = P 2 e ( 6 z )F[N^(x)][VN e e ‘\ 
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We evaluate the functional integrations in Eq. (|T(]) using a method similar to that of the lattice gauge theory. 
In the numerical calculation we discretize the neutrino path into I bins, Sxt (i = 1, 2, ...I) and in the i-th 
bin the EeD function N e (x) is given by Eq.(M). Furthermore we have replaced the functional integration 
over the EeD by a sum over K arrays, TVe k=l, 2, ...K. In Eq. © Pk is the conversion probability 
evaluated with the k'th density profile . As to PREM each point in the array N, k> which consists of 
N e (x i), ....N e (xi), ....N e (xi) is generated from PREM weighted with a Gaussian-like logarithm distribution. 
Since the deviation from PREM due to local variation is roughly 5% and the deviation is averaged per 
spherical shell with thickness of 100 Km, we take r = 5% and choose the bin sizes Sxi to be the distance the 
neutrino travels along the path of zenith angle 9 Z within a spherical shell of thickness 100 Km. So in general 
Sxi will not be equal except for 9 Z = 0. 

We note that the EeD uncertainty scale Sx differs from the one, Ip = p/ 2 f considered in 0 to characterize 
the flatness (adiabaticity) of the density profile. Both of these scales are important to the studies on the 
neutrino oscillations in matter. The effects of the l p can be taken into account in the exact numerical 
calculation, however to reduce the error caused by 5x a more precise density profile is needed. Especially 
when Sx is comparable to the neutrino oscillation length in matter, one has to be careful in estimating the 
errors for oscillation probability. 

In Fig.l.(b) we plot 5f reg and f reg as a function of the zenith angle. One can see from this figure that 
LMA suffers a larger error. For LOW the error is roughly two percent and for SMA the error is much smaller. 
So we have not shown them in the figure. Integrated over the zenith angle, it gives rise to a correction of 
20% roughly to the Ajjn for LMA, however the corrections are small for LOW and SMA. Combined with 
Fig.l.(a), we see that the errors are small so that Aj^// A and A/ffl' can be distinguished. 

To see the effects on the solar neutrino observations, we now estimate the errors in the rate of the charged 
current events during the night time. Following || we define the normalized rate of the charged current 
events as, 


[CC}(9 Z ) = Ncc/Ncc 
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where &(E V ) contains both the neutrino flux from the Boron decay and the He + P chain in the sun(l|, |f|, 
and Nqq M is the normalization factor which equals to the integral in r.h.s. taken at P N = 1. From the 2’nd 
to the 3’rd line of Eq. (0 , the integration of the differential cross section <7 with respect to the recoil electron 
kinetics T e and the scatting angle 9l has been replaced by a total charged current cross section acc of the 
neutrino on the Deuteron, since the possible uncertainty from T e , 9l can be canceled in the [CC\ as a ratio 
of Nqc to Nqq M . The E„ dependence in acc is accessed by employing a quick function from interpolation 
in |^o|. The starting point of the neutrino energy is set at E® ss Q + Eth , with Q = 1.442MeV being the 
deuterium threshold energy and Eth = 5 MeV the electron threshold energy. In Fig.2, we plot the zenith 
angle distribution of the charged current events rate in Eq. (0- One can see that the SNO charged current 
data lies in the middle between the LMA and LOW. This serves also as a check of our numerical calculation. 
Following the binning method of [||, we plot in Fig.3.(a) and Fig.3.(b) the charged current events v.s. bins 
(note for the fifth bin cos9 z ~ (0.83, 0.92) in the case of SNO), which shows that, 


LMA : [N]i < [AT] 2 < [TV ] 3 < [TV ] 4 < [TV] 5, 

LOW : [TV ] 2 > [N} 4 > [IV] 1 ~ [TV ]3 > [D\. (12) 

Quantitatively, it reads, [TV] 5 /[TV] 2 = 0.999 ~ 1 and [TV] 2 /[TV ] 3 = 0.995 ~ 1 for LMA while [TV] 5 /[TV ] 2 = 
0.982, [TV] 2 /[TV ] 3 = 1.053 for LOW, from which it might be possible to distinguish LMA from LOW. The 
calculation for SMA can be easily worked out, however for simplicity we will not repeat it here. 

Making use of Eqs. (0), ©’ and © we estimate the errors in the charged current events rate caused 
by the uncertainties in the earth electron density 


6[CC] oc f dE v $(E v )(-2X6freg)(Tcc(E v ), (13) 

Je° 

which we show in Fig.([|) by the error-bars. To avoid multi-fold integration which is a computer time con¬ 
suming, we investigate Sf reg at neutrino energies of 8 ,10,11,12 MeV and find the results almost unchanged. 
To be conservative we have used the maximal value for Sfreg- 
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We see from the figure that the errors become larger as the zenith angle increases in the case of LMA. 
Averaged over bins we have ([_/V ] 5 — < 5 [-/V] 5 )/[iV ]2 ~ 0.935 while ([N] 2 + A[] 2 )/( 3 — A[iV ] 3 ) ss 1.043. As 
indicated in the figures 13-16 of |j] that the LMA sheet in their correlation figures mainly stretched along 
the Adn direction, we study a correlation between [A r ] 5 /[A ^] 2 and [IV] 2 /[IV] 3 , which we show in Fig.5. One 
sees that the point (1, 1) for LMA is swollen into a rectangle close to the point ( 0.982,1.053 ) for LOW. In 
this figure we have not shown the error bars for LOW since they are small. 

So far the precision of PREM which we assume is 5%. Certainly errors on the zenith angle distribution 
become larger if the uncertainty in the earth electron density is bigger. Sure a modern Earth’s density model 
with higher precision will reduce the errors considered in this paper. As an example we take density model 
AK135pl||. The precision of AK135 is widely considered to be about 1 ~ 2 %, and its uncertainty scale 
is roughly 5x ~ 50 Km since the model was presented in a data table. Taking a 2% uncertainty in the 
electron density we show our results in Figure. 6 . one finds ([IV ]5 — A[iV] 5 )/([TV ] 2 + <5[IV] 2 ) « 1.033 while 
([IV ]2 + < 5 [-/V] 2 )/([-/V ]3 — 5[IV] 3 ) ~ 1.017. From Fig. 6 .(b), we see the gap between LMA and LOW enlarged. 
This makes it easier to distinguish the LMA from LOW than the prediction from PREM. 

In summary, we have estimated in this paper the errors in the zenith angle distribution of the charged 
current event rates of the solar neutrinos originated from the electron density uncertainty. Our results show 
that the corrections are not significant in the case of LOW and SMA, however, error is notable for LMA. 
Even though our estimations are given for specific parameters and qualitatively, the results of this paper 
indicate that to observe the zenith angle distribution a precise knowledge on the Earth electron density is 
necessary. 

The work is supported in part by the NSF of China under Grant No 19925523 and also supported by the 
Ministry of Science and Technology of China under Grant No NKBRSF G19990754. 
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Figure 1: Plot of the regeneration factor vs the zenith angles for neutrino energy at 11 MeV. The Earth 
matter model of PREM and the neutrino oscillation parameters in Eq.(||) have been used, (a) The solid line 
is for LMA while the dotted line for LOW. (b) the error bars corresponds to the corrections due to the 5% 
uncertainty in the matter density ( PREM ) . The fluctuation in LOW case is smaller than the LMA case, 
so we have not shown them explicitly in this figure. 
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Figure 3: Plot of the charged current event rates averaged over bins as a function of the zenith angles, (a) 
is for LMA which corresponds to the dashed line in Fig.2, (b) is for LOW corresponding to the dotted line 
in Fig.2. 
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Figure 4: Plot of the errors in the charged current event rates for LMA vs the zenith angles, (a) shows 
the error-bars attached on the dashed line of Fig.2. (b) The solid line is the same as that in Fig.3(a). And 
between the dotted lines are the errors caused by the uncertainty in the electron density. 
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